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\Q Abstract 

In this article, we prove a new functional limit theorem for the partial 
sum sequence S[ n t\ = Sl=*i corresponding to a linear sequence of the 
L. form Xi — ^2j €Z Cj£i-j with i.i.d. innovations (£i)igz and real-valued 

'"" . coefficients (cj)jgz. This weak convergence result is obtained in space 

D[0, 1] endowed with the S-topology introduced in [18], and the limit 
process is a linear fractional stable motion (LFSM). One of our result 

S provides an extension of the results of [3] to the case when the coefficients 
i (cj)jgz may not have the same sign. The proof of our result relies on 

the recent criteria for convergence in Skorokhod's Mi-topology (due to 
y—i [SDi an d a result which connects the weak S-convergence of the sum 

of two processes with the weak Mi -convergence of the two individual 
C processes. Finally, we illustrate our results using some examples and 

computer simulations. 

t-H MSC 2010 subject classification: Primary 60F17; secondary 60G52 

ON 

^ 1 Introduction 

L* The study of limit theorems for stochastic processes was initiated by Donsker in 

. £h [H] an d Prohorov in 28! in the case of processes with continuous trajectories, 

and continued by Skorokhod in his seminal article |32j . in which he introduced 
the topologies J\,M\,J^Mi on the space D[0, 1] of cadlag functions on [0,1]. 
The basic idea is that criteria for compactness for sets in D[0, 1], once trans- 
lated into criteria for tightness for probability measures on this space, become 
-via Prohorov's theorem- very powerful tools for proving functional convergence 
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of stochastic processes, in the presence of finite-dimensional convergence. Al- 
though the immediate goal of Skorokhod's work was to extend the classical 
limit theory for sums of i.i.d. random variables to functional convergence, he 
saw this as part of a bigger program in which "analysis of stochastic processes 
based on approximating them by simpler ones" plays an important role. After 
the publication of Billingsley's cornerstone monograph these ideas have been 
developed into a solid theory which has been extended in many directions, and 
nowadays has ramifications in basically every area in probability theory. 

In the present article, we study the functional convergence of the partial sum 
sequence {S n — X)™=i Xi,n > 1} corresponding to the linear sequence: 



for suitably chosen coefficients {cj)jez and an i.i.d. sequence (£;)iez- This 
problem has a very rich history and has been investigated by many authors. The 
philosophy behind these investigations is the one commonly encountered when 
dealing with a convolution between a random object (describing the shocks that 
drive the system), and a deterministic filter (the impulse-response function): if 
the filter is sufficiently smooth, then one expects that most properties of the 
input noise (in this case, the sequence (£i)igz) can be transferred to the outcome 
result (in this case, the sequence (Xj)j e z). 

It is known that for a sequence (£i)igz of i.i.d. random variables, the class 
of distributional limits of its (suitably normalized) partial sum sequence coin- 
cides with the class of stable distributions, these being the only distributions 
which possess a domain of attraction. Here we used the common terminology 
according to which a random variable £ belongs to the domain of attraction of 
a distribution G if there exist some constants a n > and b n £ K, such that 



where (£j)j>i are i.i.d. copies of £ and Z has distribution G. 

The properties of random variables in the domain of attraction of stable 
distributions depend on the value of the index a £ (0, 2] of stability, the case 
a = 2 corresponding to a normal limit distribution in ([2]). Since the constants 
a n , b n and the parameters of the distribution of Z play an important role in the 
present article, we recall below their definitions. It is important to note that 
these objects are quite different in the case a = 2, versus a £ (0, 2). 

A comprehensive unified treatment which covers simultaneously the case 
a = 2 and a £ (0,2) is given by Feller in Section XVII. 5 of [15]. We describe 
very briefly the salient points. A (non-degenerate) random variable £ is in the 
domain of attraction of a normal distribution if and only if U(x) — E(^ 2 1^ <X \) 
is slowly varying. In this case, the constants (a n ) n >i are chosen such that 




(1) 





nU{a n ) 




(3) 
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for some constant c > 0, b n — E(£)/a n and Z has a N(0,c) distribution. Of 
course, when E\£\ 2 < oo, c = Var(£). On the other hand, a random variable £ 
is in the domain of attraction of a stable law of index a £ (0, 2) if and only if 

P(|£| > 2) ~ ar"£(izi) as x — Y oo (4) 

for some slowly varying function L, and 

hm — yh-. = p and hm — — = q (5) 

x^oo P(|f| > x) Hoc P(|£| > X) W 

for some p > and 3 > with p + q = 1. In this case, b n = if a £ (0, 1), 
6„ = E(£)/a n if a € (1, 2) and 6 n = P[sin(£/a„)] if a = 1. The constants (a„)„ 
are chosen such that 

nP(|C| > a„) -»■ C (6) 

for some constant C > 0, and Z has a stable S a (a,f3,0) distribution, i.e. 
E(e iuZ ) = e^ u \u € K with 

-\u\ a a a (l-isgn(u)/?tan^) ifa^l 
— (l + isgn(w)/3| In if a = 1 



and parameters (/3, cr) given by: /3 — p ~ q and 



C^icos(^) ifa^l 
Cf ifa = l 



(7) 



Following Skorokhod's ideas, the next step in this line of investigations was 
to derive a functional analogue of pj). In the case when £ is in the domain 
of attraction of a normal distribution, this is given by Donsker's theorem if 
P(£ 2 ) < 00, respectively by Proposition 1 of pj5] when P(£ 2 ) = 00. In general, 
such a result can be deduced from Theorem 2.7 of [33] ■ More precisely, using 
the same constants a n and b n as in ([2j), one can prove that 

, [«•] 

in the space O[0, 00) of cadlag functions on [0, 00), equipped with the J\- 
topology. The limit process {Z(i)} t > is a Brownian motion if a = 2, and 
an a-stable Levy process if a € (0, 2). In both cases, Z(l) ~ 

Functional limit theorems for linear sequences of the form |TJ) with X^ez l c j 1 2 < 
co and (possibly dependent) innovations (£i)i>i with hnite variance have re- 
ceived a lot of attention in the literature, usually by treating first the short-range 
dependence case (when J2i>i E{X Xi) < 00), and then the more difficult case 
of long-range dependence (when J2i>i E(X Xi) = 00). We refer the reader to 
[TU] . [25] . [TT] for a sample of relevant references. The remaining case when £1 is 
in the domain of attraction of the normal law and has possibly infinite variance 
has been examined in the recent article 125 . 
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A totally new direction for the study of functional limit theorems for depen- 
dent sequences was initiated by Durrett and Resnick in [T3]. In the i.i.d. case, 
this method supplied a surprising new proof of Skorokhod's result ([8| based on 
the convergence in distribution of the point process N n = Yli>i^(i/n,d/a n ) to 
the underlying Poisson process N of the Levy process {Z(t)} t > (see Proposi- 
tion 3.4 of 29 ). The power of this method lies in the fact that it can be applied 
to very general dependent sequences, and supplies the convergence of a broad 
range of functionals (not only the sum). Since then, this method has been used 
by many authors in a variety of situations (e.g. [9], [17], [2], [8], [34], [4]). In 
the case of linear sequences of the form ([I]), a different method was used in [2"T] . 
by expressing the sum S n — X)"=i X% as ^ ne sum °f an i-i-d. sequence and a 
negligible term. 

From the results of [2] and [9], it can be inferred that: 

~^2Xi-A[n-}b n f ^4'AZ(-), where A = ^ Cj -, (9) 

a " »=i jez 

where denotes finite-dimensional convergence, a n , b n and Z(-) are the same 



as in (|8]), and the coefficients (cj)j £ z are such that X^jez l c il < 00 an< ^ 

\cj\ 5 < oo for some < 6 < a. (10) 

Assume that E(£) = if a 6 (1,2) and ^ has a symmetric distribution if 
a = 1, so that 6„ = 0. A natural question is to see if ([9| can be extended to the 
functional convergence 



AZ(-) (11) 



i=l 



in the space D[0, 1] equipped with a suitable topology. In the well-known article 
[3], Avram and Taqqu showed that if the series ([T]) has at least two non-zero 



coefficients Cj, then ( 11 ) does not hold in the Ji-sense. Their argument showed 
that Skorokhod's criterion for Ji-tightness is not satisfied by a linear sequence 
([I]) with finitely many non-zero coefficients. However, in the case when all 



the coefficients are non-negative, they showed that (11) holds in the Mi-sense, 
assuming for a G (1, 2), that the coefficients satisfy a technical condition (which 
holds for instance if S < 1 and (cj)i>o and (ci)i<o are monotone sequences). 



The fact that (11) cannot hold in D[0, 1] equipped with any topology for 
which the supremum is continuous can be seen more easily by considering the 
linear sequence 

Xi = £i — i e z 

whose coefficients are: Co — 1, c± — —1 and Cj = for any j E Z\{0, 1}. In this 
case, the finite-dimensional distributions of S[ nt ]/a n = (£r n t] — S,o)/ a n,t € [0,1] 
converge to 0, but snp te [ a l ]S[ nt ]/a n = max fe <„(£„ - £o)/a« converges to a 
random variable with a Frechet distribution. 
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Finding a suitable topology on O[0, 1] for which ( 11 1 holds has remained an 
open problem in the literature since the article of Avram and Taqqu in 1992. 
In the present article, we give one possible answer to this question, using the 
S'-topology introduced in [TB] (which is weaker than Ji and Mi). 

Our first result (Theorem |3.9[ ) is stated in the more general framework of [2Tj . 
in which the normalizing constants for the partial sum sequence S n are of the 
form d n a ni for regularly varying constants d n of index H—l/a with H G [1/a, 1), 
and the limiting process is a linear fractional stable motion (LFSM), which is 
an iJ-sclf similar process defined as an integral with respect to {Z(t)} t ^R. 



Our second result (Theorem 3.101 is more clearly connected with the open 
problem mentioned above, and shows that if the coefficients (cj)j^z admit a 
decomposition of the form 

c 3 = c'j - c", with Cj,d- > (12) 



for which both (c^)jgz and (c^')jgz satisfy (10) for some < S < a, S < 1, 
then (11) holds in the space O[0, 1] equipped with the S topology. By taking 
c'j = max(cj,0) and c" = max(— Cj, 0), we see that this requirement is in fact 



equivalent to condition (10); for our purposes, it is more convenient to express 



( 10 ) in this form. In particular, ( 11 ) holds (in the S-sense) for a linear sequence 
(Tj with finitely many non-zero coefficients. 

To prove this, we use the recent result of [M] which shows that for a strictly 
stationary sequence (Xi)i>i which satisfies a dependence property called asso- 
ciation (introduced in [14] ), the finite dimensional convergence of the process 
{in 1 Sl=| ■^£>* G [0, 1]} to a cadlag process Y = {Y(t)} t > is sufficient for its 
convergence in D[0, 1] equipped with M\, provided that the normalizing con- 
stants 7„ are regularly varying of index j3 <E (0, 2], and the tailsum of Y(t) is of 
order o(x~^/ 2 ). Using basic properties of the association, one observes that the 
linear sequence ([!]) is associated, if either Cj > for all j 6 Z or Cj < for all 

jez. 

The second step in our argument is to add up the partial sums S' n and S 1 " 
corresponding to the linear sequences with coefficients c^- and c", which arise 



from decomposition (12). As shown in [TS], one of the remarkable properties of 
the S'-topology is that addition is sequentially continuous. In the present article, 
we prove an analogue of this property for stochastic processes. More precisely, 
we show that for processes (X n ) n , (Y n ) n , X, Y with trajectories in D[0, 1], if the 
finite dimensional distributions of {(X n ,Y n )} n converge to those of (X,Y) and 
both sequences (X n ) n and (Y n ) n converge in distribution in B[0, 1] equipped 
with the Mi-topology, then (X n + Y n ) n converges in distribution to X + Y, in 
the space D[0, 1] equipped with the S topology. This result provides us with the 
major ingredient needed to conclude our argument, since the final dimensional 
convergence of the (suitably normalize) pair (S' n , S'^) can be deduced without 
too much effort from the results of [2"T] . 

We conclude the introduction with few words about the organization of the 
present article. In Section [2j we recall the definition and main properties of 
the S'-topology (as presented in [TB]) and we prove the new result described 
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above, regarding the convergence of a sum of two cadlag processes. The proof 
of this result relies on the fact that the ^-topology is weaker than Mi , a result 
whose proof we include here for the sake of completeness, since it was stated 
without proof in [T5] . In Section [3j we recall the definition of the LFSM and 
we derive our results regarding the convergence in D[0, 1] (equipped with the 
^-topology) of the partial sum corresponding to a linear sequence ([I]) with 
coefficients (cj)j^z with alternating signs. In Section|4j we illustrate our results 
using some examples and simulations. The Appendix contains some auxiliary 
results which are needed for the proof of the fact that the ^-topology is weaker 
than Mi. 

Acknowledgement: The authors would like to thank Pierre- Yves Gaudreau- 
Lamarre for help with some of the simulations. His undergraduate project 
entitled "Donsker theorem and its applications" was funded by the Work-Study 
Program at the University of Ottawa, during the summer of 2012. 

2 The S-topology 

In this section, we recall the definition of the .S-topology introduced in [18] and 
we prove some new properties, which are used in the present article. 

Let D[0, 1] be the space of cadlag functions x : [0, 1] — s- E, i.e. functions 
which are right-continuous on [0, 1) and have left-limits at each point t £ (0, 1]. 
The space D[0, 1] can be endowed with the uniform topology, given by the norm: 

IMU = sup \x(t)\, 
te[0,i] 

but also with the four Skorohod topologies ( J\, Mi, J%, M^) introduced in [52"] . 

The iS-topology is a sequential topology on D[0, 1], defined using the concept 
of — >s convergence, which in turn relies on the concept of weak convergence 
for functions of bounded variation. 

Let V[0, 1] be the space of functions v : [0, 1] — >• M of bounded variation, 
endowed with the topology of weak-* convergence: if (i! n ) n >i and v are elements 
in V[0, 1], we write v n — > w v if 

f(t)v n (dt)^ f f(t)v(dt) 
Jo 

for any / £ C[0, 1], where C[0, 1] is the space of continuous functions on [0, 1]. 

Note that each function v £ V[0, 1] induces a signed measure v v on [0, 1] 
defined by i/„([0,t]) = v(t) for all t £ [0,1]. Since the space of the signed 
measures on [0,1] can be identified with the dual space (C[0, 1])* of C[0, 1] 
(endowed with the sup-norm metric), the space V[0, 1] can be endowed with the 
topology of weak-* convergence of (C[0, 1])*. 

The following result is crucial for the definition of the — S>s convergence. 
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Lemma 2.1 If (v n ) n >i andv are elements in V[0, 1] such thatv n — > w v, then 
there exists a subsequence (rik)k>i & n d a countable set D C [0, 1) such that 

v nk (t) -> v{t) for all t £ [0, 1] YD. 

In addition, v n (l) — > v(l). 



The proof of Lemma 2.1 follows by Banach-Steinhauss theorem and Helly's 
compactness theorem, using the fact that any v € V[0, 1] can be written uniquely 
as v = v(0) + vi — V2 for some non-decreasing (cadlag) functions V\,v 2 with 
«i(0) = v 2 (0) = 0. We omit the details. 

We recall the definition of the S-convergence from [15] , 

Definition 2.2 Let (x n ) n >i and x be elements in D[0, 1]. We write x n — >s x 
if for every e > there exist some functions (w n ,e)n>i and w £ in V such that: 

(i) \\x n - u„, e ||oo < e for any n > 1 
(m) w n,£ — ^mi v e for any £ > 
(Hi) \\x — v e \\oo < e for any e > 0. 

The following result gives the explicit relationship between the sequence 
(%n)n>i and its limit x. Its proof is based on Lemma 2.1 and a diagonal argu- 
ment. We omit the details. 



Lemma 2.3 Let (x n ) n >i and x be elements otD[0, 1]. If x n — >$ x, then there 
exists a subsequence (nk)k>\ and countable set D C [0, 1) such that 

x 7lk (t) x(t) for all t e [0, 1]\D. 

In addition, x n (l) — > x(l). 

The space D[0, 1] endowed with the — >s convergence is of type £, i.e. 

(i) if x n = x, n > 1 is a constant sequence, then x n — >s x '- 

(ii) if x n — >s x then x 7lk — >s for any subsequence (nk)k- 

Therefore, we can define the sequential topology induced by the — >s con- 
vergence. A set F C D[0, 1] is S'-closed if for any (x n ) n C F with x n — ^5 x, we 
have x € F. The 5-topology is the collection of all sets G c D[0, 1] such that 
G c is S'-closed. 

We denote by — >g the convergence in the S-topology: x n — >s % if for any 
S-open set G, there exists N > 1 such that x n £ G for all n > N. By the 
Kantorovich-Kisynski criterion, x n — >s % if and only if for any subsequence 
(^fe)fc there exists a further subsequence (rife,); such that x rik[ — >s x. 

Remark 2.4 Note that, if (x n ) n >i and x are elements in D[0, 1] such that: 

(i) (x n ) n >i is S-relatively compact; and 

(ii) for any subsequence (nk)k>i with x nk — >s V we have y = x, 
then x n %■ 



7 



The most important property of the S-topology is the characterization of its 
relatively compact sets. This is expressed in terms of the number of upcrossings, 
or the number of oscillations, whose definitions we recall below. 

For any real numbers a < b, let N a ' b (x) be the number of upcrossings of the 
interval [a, b] by the function x £ D[0, 1], i.e. the largest integer N > 1 for which 
there exist some points 

< h < t 2 < t 3 < t 4 < . . . < *2W-1 < t 2N < 1 (13) 

such that Xt 2k _ 1 < a and x t2k > b for all k = 1, . . . , N. 

For any rj > 0, let N v (x) be the number of n-oscillations of the function 
x £ D[0, 1], i.e. the largest integer N > 1 for which there exist some points 
(ti)i<i<N satisfying (13), such that 

\x(t 2k ) - x{t 2k -i)\ >V for all k = 1,...,N. 

Note that for any function x G D[0, 1] and for any a < b, 

N a > b (x) < N b _ a (x) < oo. 

Lemma 2.5 (Lemma 2.7 of [18j ) A set A c D[0, 1] is S -relatively compact 
if and only if it satisfies the following two conditions: 

(i) sup Hxlloo < oo 

(ii) sup N a ' h (x) < oo for all a,b el,s < b. 
Condition (ii) can be replaced by: 



(ii)' sup N v (x) < oo for all n > 0. 

The relationship between the S'-topology and the Mi-topology plays an im- 
portant role in the present article. We refer the reader to the original Skorohod's 
article 32 for the definition of the Mi topology, as well as Chapter 12 of [55] 
for a comprehensive account. 

The Mi -convergence can be described using the oscillation function: 

w{x,S):= sup H(x(ti),x(t 2 ),x(t 3 j), 

OV(t 2 -<5)<ti<t 2 <t3<lA(t 2 +i5) 

for any x £ D[0, 1] and <5 > 0, where H (a, b, c) is the distance between b and the 
interval with endpoints a and c: 

H(a, b,c) — (oAc-aAcAi))V (aVcVd-aVc). 
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Theorem 2.6 (2.4.1 of [32]) Let (x n ) n >l andx be arbitrary elements in O[0, 1]. 
The following two conditions are necessary and sufficient for x n — >m x x: 

(a) x n (t) x(t) for any t € Q, for a dense set Q C [0, 1] containing and 1; 
(6) lim limsup w(x n , S) = 

For the necessity part, Q could be the set of continuity points of x, together with 
and 1. 

The following result was stated without proof in [18 . A short proof can be 
given using Skorohod's criterion 2.2.11 (page 267 of 32J) for the Afi-convergence, 
expressed in terms of the number of upcrossings. This proof has a clear disad- 
vantage: it refers to an equivalent definition of the Mi -convergence, which was 
not proved in Skorokhod's paper. In the present article, we give a new proof. 

Theorem 2.7 The S -topology is weaker than the M\-topology (and hence, weaker 
than the J\-topology). Consequently, a set A C D[0, 1] which is Mi-relatively 
compact is also S -relatively compact. 

Proof: Let (ai n )„>i and x be elements in D[0, 1] such that x n — >M\ x - We will 



prove that x n — >s %■ For this, we apply Remark 2.4 Suppose that (rik)k>i is 
an arbitrary subsequence such that x„ k — >s y. By Lemma [2. 3[ there exists a 
subsequence (n^,); such that x nk (t) — > y(t) for any t G Q\ for some dense set 
Qi C [0,1] which contains 1 (whose complement is countable). On the other 



hand, by Theorem 2.6 x„ k (t) — > x(t) for any t G Q 2 for some dense subset 
Q2 C [0,1] which contains and 1 (whose complement is countable). Hence 
y{t) = x(t) for any t £ Qi D Q2, and y = x. 

"it 



2.5 



It remains to prove that (x n ) n >i is S'-relatively compact. By Lemma 
suffices to prove that: 

sup HaCnlloo < 00 (14) 
n>l 

sup N v (x n ) < 00 for all i] > 0. (15) 

n>l 

Since x n — >m x % and the supremum is Mi -continuous, ||a; ra ||oo ~~ * Halloo- 



Relation ( 14 1 follows 



The remaining part of the proof is dedicated to ( 15 ). Let rj > be arbitrary. 



Let < e < 77/2. By Theorem 2.6 there exist some 5 > and an integer 
no > 1 such that 

w(x n ,S) < e for all n > n . (16) 

Since the set Q of continuity points of x is dense in [0, 1], we can find some 
points = to < t\ < . . . < tjtf = 1 in Q such that for each j = 0, 1, . . . , M — 1, 

t j+1 - tj < S. 

By Theorem |2.6[ there exists an integer n\ > hq such that for any n>n\ 
\x n {tj) - x{tj)\ < e, for any j = 0,1,..., M. (17) 
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Fix an integer n > n\. Suppose that there exist some points 

< Sl < S 2 < S 3 < S 4 < . . . < S 2 N~1 < S2N < 1 

such that 

\x n {s 2 k) - x n (s 2 k-i)\ > V> for all fe = 1,2, ...,iV. (18) 



The proof of (|15j) will be complete once we estimate the number TV by a constant 
independent of n. 

Relation ( 18 ) says that the function x n has N ^-oscillations in the interval 
[0,1]. These oscillations can be divided into two (disjoint) groups. The first 
group (Group 1) contains the oscillations for which the corresponding interval 
[s2fc-i, S2k) contains at least one point t' y Since the number of points tj is M , 

the number of oscillations in Group 1 is at most M. (19) 

In the second group (Group 2), we have those oscillations for which the 
corresponding interval [s2k-i, s 2 k) contains no point tj, i.e. 

tj < s 2 k-i < S2k < tj+i for some j = 0, 1, . . . , M - 1. (20) 



We now use Lemma A. 2 (Appendix A). Note that 

/3 n := sup H(x n (u),x n (v),x n (w)) < w(x n ,S) < e, 

tj <U<V<W<tj+i 

and hence, 

i r, , i\ ^ ^\ x n{tj + l) x n{tj)\ "i" Pn sup + E 
^ri{ x n, + l I S a < 

where N v (x n , [tj, tj+i]) is the number of ?y-oscillations of x n in the interval 
[tj,t j+1 ] and 

K sup = Sup Halloo < +00. 
n>l 

Since there are M intervals of the form [tj,tj+i], we conclude that 

AK SU + e 

the number of oscillations in Group 2 is at most M ■ — (21) 

r\ — e 



Using (19) and (21 1, we conclude that 



AK sup + e\ n/ A K su P + V 



N < M 1 + — = M 



r] — e J t] — e 



which does not depend on n. This concludes the proof of (15). □ 

We now provide an example of a sequence (a; n ) n >i in D[0, 1] which is S- 
convergent, but does not converge in the Mi topology. 
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Example 2.8 Let x — and 



%n(t) = l[l/2-l/n,l](i) - l[l/2 + l/n,l](i) 



1 

otherwise 



Then x n — >s x. To see this, we take v Ut£ = x n . Then i?„ j£ — ^ u e = since 
forfuiy/€C[0,l], 

'/(*)*.«) -/(i + i)-,(i-i)-,o. 



The fact that (x n ) n >i cannot converge in M\ follows by Theorem 2.6 since if 
h < I - i < h < I + i < * 3 , then fl"(K n (*i),a: n (t2),Xn(*3)) = 1- 

We now consider random elements in D[0, 1] endowed with the S-topology. 
Recall that the Borel er-field generated by the ^-topology coincides with the 
Kolmogorov's er-field V generated by the projections nt,t € [0,1] (see [T5]). 
Therefore, a random element in D[0, 1] is a random variable X : fi — > D[0, 1] 
which is measurable with respect to P. Its law is a probability measure on 
),1],2?). 



The following result is an immediate consequence of Theorem 2.7 



Corollary 2.9 Let (P a ) a <£A be a family of probability measures on D[0, 1]. If 
(P Q ) Q gA is uniformly M\-tight, then (P Q )aeA is also uniformly S-tight. 



By Lemma |2.5[ one can give some criteria for the S'-tightness of a family of 
probability measures on D[0, 1]. More precisely, we have the following result. 

Proposition 2.10 (Proposition 3.1 of [18j) Let (X a ) a£ \ be a family of ran- 
dom elements in D[0, 1]. The family (X a ) ae \ is uniformly S-tight if and only 
if it satisfies the following two conditions: 

(i) {|j^a||oo}aeA is uniformly S-tight 

(ii) {N a ' b (X a )} ae \ is uniformly S-tight, for all a, b € M, a < b. 
Condition (ii) can be replaced by: 

(ii)' {A^(X q )} q6 a is uniformly S-tight, for all 77 > 0. 

Despite the fact that D[0, 1] equipped with the ^-topology is not a metric 
space, the Direct Half of Prohorov's Theorem still holds, i.e. a uniformly S-tight 
family probability measures on O[0, 1] is relatively compact with respect to the 
S-weak convergence (defined using the convergence of integrals of ^-continuous 
functions on D[0, 1]). This follows by a strong form of Skorohod's Representation 
Theorem (Theorem 1.1 of [12]), using the fact that the S-topology possesses a 
countable family of 5*-continuous functions which separate the points in D[0, 1]. 

The Converse Part of Prohorov's Theorem also holds, but for this one needs 
to consider a stronger form of convergence in distribution (denoted by — >x>), 
given by the following definition (which was originally introduced in [20]). 
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Definition 2.11 Let (X n ) n >i be random elements in D[0, 1]. We say that 
(X n ) n >i is ^^p-convergent in distribution if for any subsequence (n k ) k 
there exists a further subsequence (n kl )i such that 

(X nh )i admits a strong a.s. Skorohod representation, 

i.e. there exist some random variables (Yz)z>i and Y defined on ([0, 1], B[0, 1], A) 
with values in D[0, 1] such that: 

(i) Yi has the same distribution as X„ k[ for any I > 1; 

(ii) y,(w) -^s Y(u) for any u £ [0, 1]; 

(iii) for any e > there exists an 5-compact set K C B[0, 1] such that 

P(Yi G K for all I > 1) > 1 - e. 

Remark 2.12 Let g : D[0, 1] — > K be an S-continuous bounded function. As 
a consequence of the previous definition, we obtain that E[g{X n )\ — > E[g(Y)], 
since for any subsequence (n k ) k there exists a further subsequence (n kl )i for 
which E[g[X nu )] — ¥ E[g(Y)\. This proves that (//„)„ converges weakly to \i 
(with respect to S), where [i n is the law of X n and /x is the law of Y. If X is a 
random element in D[0, 1] with law /i, we write J n — >%> X. 

The following two results are needed for the proof of Theorem |2.15| below. 
We recall them for the sake of completeness. 

Theorem 2.13 (Theorem 3.4 of [18j) Let {X a ) ae tv be a family of random 
elements in O[0, 1]. Then (X a ) a£ \ is uniformly S -tight if and only if it is 
relatively compact with respect to -^t>- 

Theorem 2.14 (Theorem 3.5 of |18j) Let (X n ) n >i and X be random ele- 
ments in D[0, 1] such that: 

(i) (X n (ti),...,X n (t k )) — >v (-X"(ii),...,X(i fc )) for any t x < . . . < t k in Q, 
and k > 1, where Q is a dense set in [0, 1] which contains 1; 

(ii) (l„)„>i is relatively compact with respect to — >x>- 

Then X n — "~^x> X. 

We conclude this section with a new result which shows that the sum of 
two processes converges in distribution (in the sense of — provided that the 
two processes converge in Mi and the finite-dimensional distributions converge. 
This result will be used in Section [3l 

Theorem 2.15 Let (X n ) n >i, (Y n ) n >i, X andY be random elements mO[0, 1] 
such that: 

(i) X n — >x> X in D[0, 1] equipped with Mi; 

(ii) Y n — >x> Y in -D[0, 1] equipped with Mi; 

(in) (X n (ti),Y n (ti), . . .,X n (t k ),Y n (t k )) — (X(ti),Y(ti), . . .,X(t k ),Y(t k )) 
for any < t± < . . . < t k < 1 and k > 1. 
Then X n + Y n X + Y. 
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Proof: First, let us observe that both sequences (X n ) n >i and (Y n ) n >i arc 
uniformly Mi -tight. This follows by Le Cam's theorem (Theorem 8 in Appendix 
III of [5]), since a single probability measure on D[0, 1] is Mi-tight. 

By Corollary 2.9 both sequences (X n ) n >i and (Y n ) n >x are uniformly S- 
tight, and hence, they satisfy conditions (i) and (ii)' of Proposition 2.10 Since 
for any functions x,y € D[0, 1], ||x + y||oo < ||a;||oo + IMloo and 

N v (x + y) < N v/2 (x) + N v/2 (y) for any 77 > 0, 

it follows that the sequence (X n + Y n ) n >i also satisfies conditions (i) and (ii)' 



of Propos ition 
Theorem |2. 13 



2.10 



(X n 



Hence, the sequence (X n + Y n ) n >i is uniformly 5-tight. By 
+ Y n ) n >i is relatively compact with respect to — -)*z>. 



Finally, by our hypothesis (in) and the Continuous Mapping Theorem, 

(X n (ti) + Y n (h), X n (t k ) + Y n (t k )) (X(h)+Y(ti), . . . , X(t k ) + Y(t k )) 

for any Q < t± < . . . < tf~ < 1 and k > 1. The result follows by applying 
Theorem 2.14 to the sequence (X n +Y n ) n >\. □ 



3 Our Results 

In this section, we derive some new results regarding the convergence in distri- 
bution (in the 5-topology) of the partial sum of the linear sequence . 

Let (X,-),- e g be the linear sequence given by ([!]), where (£i)iez are i.i.d. 
random variables in the domain of attraction of a stable law of index a € (0, 2], 



and the coefficients (cj)j e z satisfy (10) 



As in [3], we assume that E(£i) = if a G (1,2] and £1 has a symmetric 
distribution if a = 1. Note that the series appearing in the right-hand side of 
([I]) converges a.s. This follows by applying Theorem 12.10.4 of [22] (when (10) 



holds with S > 1), or Theorem 12.11.2 of [22] (when (10) holds with 8 < 1) 



For any n > 1, we consider the partial sum processes: 

[nt] [nt] 

Z n (t) = —Yji,t>0 and S n (t) = yX l ,t>0. 



i—i 



Due to our assumptions, relation ^ holds with b n — 0, i.e. 

Z n (-) -^V Z(-) 

in D[0, 00) equipped with the Ji-topology. The constants (a n ) n and the distri- 
bution of Z(l) ~ Z have been specified in the introduction. 

Remark 3.1 For any a € (0,2], Z has a strictly a-stable distribution. When 
a = 2, {Z(i)} t >o is a Brownian motion with variance c. In this case, Z(t) —Z(s) 
has a N(0, c(t — s)) distribution, for any s < t. 
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When a <E (0,2), {Z(t)} t >o is an a-stable Levy motion, in the sense of 
Definition 3.1.3 of |31j (except that in this definition, it is required that a = 1). 
To see this, note that for any s < t, by property 1.2.3 (page 11 of [31] ) 

2(o-^)~^-.)~(.-.)^~{*;w;::p^ ) 



where /ii = —^@-(t— s) 1 /" ln(t— s). If a = 1, we assume that p — q and fii = 0. 

We follow very closely the approach of [21] . We first observe that a~ 1 S'„(t) 
can be expressed as an integral with respect to the process Z„ = {Z n (t)} t >o. 
Recall that the integral of a deterministic function / : K — > K with respect to 
Z„ is defined by: 



/OO i 
f(u)dZ n (u) : =-J2f(j/n)^, 



provided that the sum converges a.s. One sufficient condition for this is: 



|/(j/n)| b < oo for some < b < 



(22) 



(23) 



jez 



In our case, using ([I]) and Fubini's theorem, we obtain that 

[nt] I [nt]-j \ 

Sn(t) = EE " E E C fe Kj> 

i=i jgz jez \fc=i-j / 



and hence 



where 



—S n (t) = — ^2fn (t, - J = / f n (t,u)dZ n (u), 



[ni]-[nu] 
fn(t,u) = ^ c k, UE 
k—l—\nu] 



Remark 3.2 Note that condition (23) satisfied for / = f n (t,-). This follows 
by taking b = 8 since 



E 

jez 



[«t]-i 

E * 

fc=i-j 



[nt]— j [ni] — 

jezk=i-j keij=i-k 

= C n , tiS [nt]J2\ c k\ S < oo 

fcez 



where C„ )ti <5 = 1 if S < 1 and C n ,t,$ = [nif 1 if 5 > 1 
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In [2T], it is shown that the finite dimensional distributions of the process 
S n /{d n a n ) converge to those of a LFSM, for a sequence (d n ) n of suitable con- 
stants. We will need this result below. To introduce the LFSM in the case 
a G (0, 2), we have to recall the definition of an a-stable random measure. 

Definition 3.3 Let m be a positive measure on (R, B) and Bq = {A G B; m(A) < 
oo}. Let a G (0,2). A collection M = {M(A);A G B } of random variables 
defined on a probability space (CI, J- ', P) is called an a-stable random measure 
on (R, B) with control measure m and skewness intensity f3 € [—1,1] if: 

(i) for any disjoint sets A\, . . . , Af. G Bq, M(A\), . . . , M(Ak) are independent; 

(ii) for any disjoint sets (A n ) n >i C B with U n >iA n G B , M(U n >iA n ) = 
En>i M (A n ) a.s.; 

(iii) for any A G 6 , M(A) - ^(m^) 1 /", /3, 0). 

The existence of M is shown in Section 3.3 of [3T]. The stochastic integral 

1(f) = I f(u)M(du) 



is defined in Section 3.4 of [31] for any measurable function / : R — > R such that 
/ K |/(u)| a m(<iu) < oo if a 7^ 1, and L In \ f(u)\\m(du) < oo if a = 1. By 
Property 3.2.2 (page 117 of [3T]), 

I(f)~S a {v f ,l3f,ii f ) 

where cr^ = / R |/(u)| Q m(dw), ^ = Jr/( u ) ln \f(u)\m(du) if a = 1, 

n - f A a P j\f(u)\ a sgnf(u)m(du) 

H f = if a^l and f3 f = — . (24 

J |/(u)| a m(dtt) 

In what follows, we assume that m — o~ a \ where A is the Lebesgue measure 
on R and a is given by 0. Then 

M((s,t}) ~ S a (cr(t- s) 1/a ,/3,0) for any s < t. 

Setting Zm(£) '■— M([Q,t]),t £ R, we see that {2'M(i)}t>o has the same 
finite-dimensional distributions as {Z(t)} t >o- For this reason, we say that Zm 
is an extension of Z to R. We denote Z M simply by Z and Lf(u)M(du) by 
L f(u)dZ(u). We are now ready to give the definition of the LFSM. 

Definition 3.4 (Definition 7.4.1 of [31]) The linear fractional stable motion 
(LFSM) is the stochastic process {A. a ^^ a ^{t)}t^R given by: 

/oo 
f a ,H,a,b(t,U)dZ(u), 
oo 

where a € (0,2], H G (0,1), H f 1/a, a G R, 6 G R with |a| + |6| > and 

/a,H,a,6(t, u) = o{(t - n)^ 1/Q - (-«)^ 1/Q } + 6{(u - t)^ 1/Q - («)^ 1/a }- 
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Remark 3.5 (i) The fractional Brownian motion (FBM) of Hurst index H £ 
(0, 1) is a zero-mean Gaussian process {Bh (*)}t>o with E[Bh (t)-Bff(s)] = (t 2H + 
g 2H — |t — s\ 2H )/2. This process admits the "moving average" representation: 
(see Proposition 7.2.6 of [3Tj ) 

/oo 
{(t- u )f- 1/2 -(-^- 1/2 }ds( w ) 
-00 

where {B{t)}teR is a two-sided standard Brownian motion and 

C H - 



-J 1/2 



(l + ^- 1 / 2 _ u H-l/ 2 ] 2 du+ _^| (25) 

Therefore, the process h.2 t H,c H $ is a FBM. Note that the sample paths of the 
FBM are 7-Holder continuous, for any 7 <G (Q,H). 

(ii) It is convenient to have a unified notation which covers also the case 
H = l/a. Therefore, if H = l/a we let f a ,H,a.b{t,u) = ( a ~ b)l(o,t](u) and 

/oo 
f a ,H,a,b( t > u )^' Z ( u ) = ( a ~ b)Z(t). 
-00 

In this case, A a ^n,a,b has a cadlag modification (see e.g. Theorem 5.4 of [50]). 

(hi) Suppose that a <G (1,2) and H > l/a. By Proposition 7.4.2 of [3T] . 
Aq. h a ,6 is an -ff-sssi process, i.e. it is iJ-self similar and has stationary incre- 
ments. Note that 

A- a ,ff, a ,b(t) = 1(f) ~ S a (a f ,f3 f ,fi f ) 

with / = f a ,H,a,b and parameters (<7/,/3/,/i/) as above. By (24| and our as- 
sumption that P = if a = 1, it follows that A ai £r jaj b(i) has a strictly a-stable 
distribution. By Theorem 12.4.1 of [31], A a jj a b has a continuous version. 

The following recent result of [24] lies at the origin of our investigations. We 
recall it for the sake of completeness. 

Theorem 3.6 (Theorem 1 of |24|) Let (Xi)i>i be a strictly stationary se- 
quence of associated random variables and S n (t) — X)!=*l^i,£ ^ [0, 1] . Let 
(7n)n>i be a non- decreasing sequence of constants such that 7„ — > 00 and 

7n ~ n^L(n) 

for some /3 € (0, 2] and a slowly varying function L. Let (b n ) n >\ be a sequence 
of real numbers such that K := swp{\bk — b n \; n > 1, n < k < 2n} < 00. Assume 
that either K — or liminf„ j n /n > 0. Let 

Y n {t) = — (S n (t)-[nt]b n ) te[0,l]. 

IfY n (-) Y(-) where Y = {Y(t)} te mx\ is a cadlag process with 

lim ^ /2 P(|F(i)| > x) = /or a;/ t e [0, 1], (26) 
then F„0) Y(-) 2 n (D[0,l],Mi). 
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Based on Theorem 3.6 we can derive our first result, which is an extension 
of Theorem 5.1 of |21) from the convergence of finite-dimensional distributions 
to the convergence in distribution in D[0, 1], in the Mi-sense (in the case of 
non-negative coefficients). For this, we need to observe that 

a„~n 1 / Q J L 1 (n) (27) 

for a slowly varying function L\. 

Theorem 3.7 Let a 6 (1, 2]. Suppose that Cj > for all j E Z and 
1 n 1 ° I d 

j—0 j — — n x 1 1 

for some a, b E M. with \a\ + |6| > 7 where 

d n = n H -^ a L 2 (n) (28) 

for some 1/a < if < 1 and a slowly varying function L 2 - ie^ (fl„) n 6e chosen 
to satisfy or |(|) depending on whether a = 2 or a G (0, 2). T/ien 

-^Sn(-) — Aa,i? i0l 6(.) in (B[0,l],Mx). 

Proof: To simplify the notation we write A instead of A a u.a ,b- By Remark 
|3.5[ we may assume that A is a cadlag process. We apply Theorem |3.6| with 
b n = 0, 7„ = d n a n and j3 = 1/if. By Theorem 5.1 of [2"T] . 



d n a 



Due to (27l and (28), c„ ~ n H L(n) for the slowly varying function L := LiL 2 - 
If a = 2, A(i) has a normal distribution, whereas if a £ (0,2), A(t) has an 
a-stable distribution. In both cases, (26) holds. 

Finally, we prove that (Xj)j 6 z is an associated sequence. By definition, 
X, L = limjv^oo X> N > a.s. where x[ N) = Y^=-N c i^i-r % property (P5) of 
[2], it suffices to show that (-^Q- )i£Z is associated for any N > 1. Note that 
A 4 (Ar) = /(&-Ar,...,&+iv), where the function / : R 2Ar+1 -> R is defined by 

f(xi-N, ■ ■ ■ ,%i+N) — Y^,jL-N c j x i~j- Since Cj > for all j e Z, the function / 
is coordinate-wise non-decreasing. Let i\ < . . . < ik be a finite set of indices in 
Z. Let I = {ix — N, . . . , ik + N}. Then for each j £ {1, . . . , k}, we can say that 
Xj. — fj(&,i G i) for some coordinate- wise non-decreasing function fj. Since 
(£i)iei are associated (Theorem 2.1 of [H]), by property (P4) of [H], it follows 
that {X;\ . . . , X.- } are associated. □ 



Remark 3.8 (i) Under the assumptions of Theorem 3.7 if if > 1/a then 

J^— SnG) Aa.H.a.lO) in (D[0,1],Z7) 
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since {Aa,H,a,b(t)}te[o,i] * s continuous. Here U denotes the uniform topology. 
When H = 1/a, Theorem 3.7 can be viewed as a variant of Theorem 2 of [3]; 
in this case, the limit is (a — b)Z(-). 

(ii) In the case H < 1/a, Theorem 5.2 of [H] gives the finite-dimensional 
convergence of (a^S,^-) — AZ n (-))/d n where A = J2j e %Cj. However, this 
process cannot be identified with the partial sum process corresponding to a 
sequence of associated random variables. Therefore, the case H < 1/a cannot 
be treated by the methods of the present article. 

The following theorem is the main result of this article. 

Theorem 3.9 Let a E (1,2]. Suppose that Cj = c'j — c" for some c'j,c" > 
such that dj + c" — 0(d\j\/\j\) as \j\ — > oo and 

(29) 

for some a',a",b',b" e K with \a'\ + \b'\ > and \a"\ + \b"\ > 0, where d n 
satisfies (28) for some 1/a < H < 1. Let a = a' — a" and b = b' — b" . Let (a n ) n 
be chosen to satisfy or depending on whether a — 2 or a <E (0, 2). Then 

j >Sn(0 A a) fl- iCti {,(-), 



where — >t> denotes the convergence in distribution specified by Definition 
In particular, 

t^S„(-)^A„, w (.) in (D[0,l],5). 

Proof: For any i G Z, we define X| = ^2j € z c j^-3 an ^ ^i' = Sjez^^-i* 
Then S B (t) = S;(t) - S«(t), where S;(t) = Effl and S;'(t) = Ed XV. 
By Theorem [3~7] 

1 



l 



S^(.) — > D A'(-) in (D[0,1],M!) 



S^(.)— ^A"(-) in 00(0,1], Mi), 

where A'(-) = A a>H , a /,y(-) and A"(-) = A a , H ,a",b"{-)- 

By the linearity property of the stable integrals (page 117 of [3"T]). 

/•oo 

A'(t)-A"(i)= / [f a ,H,a',b'(t,u)-f at H,a",b"(t,u)]dZ(u) = A atH ,a,b(t) a.s. 



The conclusion will follow from Theorem 2.15 once we prove that for any 
< t x < . . . < t k < 1, 

-^—{S' n {h), S'M, S' n (t k ), S^tk)) -^v (A'(h), A'%), A'(t fe ), A" (*,)). 

(30) 
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Note that each of the components of the random vector on the left-hand side 
of (30) is a stochastic integral with respect to Z n , in the sense of (22 1. More 
precisely, 

1 



d n (x ri 



S' n (t)= / f' n (t,u)dZ n (u), 



fZ(t,u)dZ n (u) 



with 



[nt] — [nu] [nt] — [nu] 

fnM = E C 'fe and /«(*.«) = E d k- 



fc=l- 



d n 



fe=l— [nu] 



The components of the random vector on the right-hand side of ( 30 ) are stochas- 
tic integrals with respect to {Z(t)} te ^. To prove (30), we apply Corollary 3.3 
of [H] to the functions: 

fl(u) = f^t 1 ,u)J^u) = f , :(t l ,u),...J^- 1 (u) = f n (t k ,u), 

ti k (u) = fZ(t k ,u), 

f l {u) = J a M,a\b'{tl,u) 1 f 2 (u) = fa,H,a",b"(tl,U), . . . , /^^(u) = f a ,H,a',b'(tk,U), 

f k (u) = f a ,H,a",b"(tk,u). 

As in the proof of Theorem 5.1 of [51], one can show that for any t £ [0, 1] 
fixed, the functions f„(t, •) and f"(t, ■) satisfy the following conditions: 
(Bl) fn(t,u) f a ,H,a',b'( t > u ) and /n(*> u ) ^ fa,H,a",b"(t,u) a.e.(ii). Here, 
— ^4 denotes the continuous convergence, i.e. f n (u) — — > f(u) at u — uq if 
fn(u n ) -> /(ii ) whenever u„ -> u . 
(B2) For any T > there exists j3 > a such that 



SU P/ l/4(^ u )l /3 Pn(dw) < oo and sup/ |/£(t, u)fp n (du) 

n>lJ\ u \<T ri^lJlul^T 

where p n is the measure on R defined by /5 n ([0, it]) = [nu]/n. 
(B3) There exists e > such that 



< oo, 



lim limsup 



(\f n (t,u)\ a - e + \&(t,uT +s )dp n (u) = 



\u\>T 



lim limsup 



\u\>T 



(X(t,uT-e + \rt(t,uT + e)dPn(u) = 



Therefore, the functions /*, . . . , /„ , f 1 , . . . , f 2k satisfy the conditions (Al)' , 
(A2)' and (A3)' of Corollary 3.3 of [ft]. Relation ([30]) follows. □ 

Our final result gives a set of conditions on the coefficients (cj)jez, for which 
we can give an answer to the open problem mentioned in the introduction. 
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Theorem 3.10 Let a G (0, 2]. Suppose that Cj = c'j — c" where c'-, c" > 0, 
^( C ;-) 5 < oo and ^(c^) 5 < oo 



.ie- 



for some < 5 < a, 5 < I. Let A = X^ez c j an d ( a n)n be chosen to satisfy 
or (ml depending on whether a — 2 or a £ (0, 2). Then 



Sn(-) 



AZ(-) in (D[0,l],5), 



where {Z(t)} t >o is the process described in Remark 



3.1 



>x> denotes the 



convergence in distribution specified by Definition \2.11 In particular, 
1 



-S n {-) — >t> AZ{-) in 



MS). 



Proof: The argument is similar to the one used in the proof of Theorem |3.9| 
The idea is to first prove the result for non-negative coefficients (cj)jgz (similarly 
to Theorem 3.7 1, and then apply Theorem 2.15 To prove the finite dimensional 
convergence of the pair (S' n /a n , S'^/a n ), we use again Corollary 3.3 of [3T] (as 
in the proof of Theorem 5.3 of |21|). □ 



4 Examples and Simulations 

In this section we consider several examples to illustrate the results presented in 
Section [3] We let £, (£i)iez be i.i.d. random variables with a symmetric Pareto 
distribution with parameter a G (0, 2], i.e. £ has density 

f( x ) = ^a^~ Q ~ ll (i,oo)(a;) + ^"(-a:) _a_1 l(-oo 

If a > 1, E(£) — 0. Note that |£| has a Pareto distribution with parameter a, 
i.e. |£| has density g(x) = ax~ a ~ 1 l^i^ 00 ){x). 

When a < 2, we choose a n such that nP(|£| > a n ) — > 1. A suitable choice is 

a n = M{x > 0;P(|f| < x) > 1 - n" 1 } = - n' 1 ) = n 1 ' 01 

where G(x) = P(|£| < x) = 1 - x~ a and G~ x {y) = (1 -y)- 1 ^^ G [0,1]. In 
this case, {.Z(i)} t >o is an a-stable Levy motion with 2(1) ~ ^(cr, 0, 0) and the 
scale parameter a is given by Q with (7=1. 

When a = 2, we choose a„ as the largest root of the equation x 2 = 2nlnx, 
since the function U(x) — E(\^\ 2 l^\<: x y) = 21ncc is slowly varying, and hence £ 
is in the domain of attraction of the normal law. (However, note that -E(£ 2 ) = 
oo.) In this case, {Z(t)} t >o is a standard Brownian motion. More precisely, we 
used the formula a n = exp{— I/ 71 )} where W—\ is the second branch of 

the Lambert W function (see e.g. [5]). 
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For the simulations, we used 



& = F- 1 ^) = -(2C/ J )- 1/Q l [0 ,i/2](^) + [2(1 - [/i)r 1/a l(i/ 2 ,i](E/i) 

where (f/i)igz are i.i.d. random variables with a uniform distribution on [0, 1], 
and -F -1 is the generalized inverse of the c.d.f. F of £: 

2- 1 (-x)- Q ifx<-l 
F{x) = { 2" 1 if x G (-1,1] 

l-2- 1 x- Q ifx>l 

We used the truncated series = 5Z|j|<jv c j£i-j as an approximation for Xj, 

and the corresponding partial sum sequence S„ (■) = Si=i -^/^ as an approxi- 
mation for S n (-), if JV is large. We considered n — 1000 and N — 50. Therefore, 
a n = 95.4883 when a = 2. 

Remark 4.1 (i) Usually, to illustrate a functional limit theorem by the conver- 
gence of plots, one uses the automatic scaling done by the computer, as explained 
in Section 1.2.4 of [3S]. For instance, for illustrating Donsker's theorem for the 
partial sum sequence S„(-) = X)l=i£i associated with i.i.d. random variables 
(£i)i>i with mean and variance 1, one would like to plot the step function de- 
termined by the points (k/n, Sk/y/n) for k = 1, . . . ,n, where S n — S n (l). To do 
this and guarantee that the plot fits in the available space, the computer auto- 
matically shifts the values Tk, n = Sk/y/n by minfc< n Tfe in and then scales them 
by a factor equal to the range r = maxj-< n Tfc >n — min/j<„ Tfc j7l . Since the map 
plot : D[0, 1] — > D[0, 1] defined by plot(x) = (x — inf tG [ 0i i] x(t))/(sup tg j 01 i x{t) — 
inf tg r 01 i x(t)) is Ji-continuous, by the continuous mapping theorem, 

plot plot(Z(-)) in (D[0, 1], Ji), (31) 



where {Z(t)} t ^ i] is a standard Brownian motion. Therefore, the convergence 
of the plots is preserved, despite the automatic shift-and-scale operation. 

(ii) In our case, we had to avoid this automatic shift-and-scale operation since 
the map plot is not ^-continuous. However, the map plot a j, : D[0, 1] — > D[0, 1] 
defined by plot a j,(x) = (x — a)/(b — a) is 5-continuous for any a,t € I with 
a < b. Therefore, for the simulations in Example |4.4| below we chose ourselves 
the range for the values Tj? n = /(d n a n ) with k = 1, . . . , n, where S% = (1) 
(for n and N fixed). More precisely, we forced a shift-and-scale of these values 
in a more "realistic" range [ir m i n , x max ]. This range was obtained by simulating 
(TJ^ n )k<n a large number M of times, recording the values T m ; n = minfc< n T^ n 
and T max = maxi<„ n each time, and taking x m \ n (respectively x max ) as 
the 10%-quantile of the M values T m ; n (respectively the 90%-quantile of the 
M values T max ). We considered M = 75. The same procedure was employed 



in Example 4.5 for = S%? /a n , using the 15%-quantile, respectively the 

85%-quantile. 
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Remark 4.2 (i) The procedure that we explained in Remark 4.1 (i) provides 
some reasonable justification for the illustration of Donsker's theorem using the 
plots. However, there is a small problem with this justification, since relation 



(31 ) gives a convergence in distribution, and not an a.s. convergence. In reality, 



what we illustrate with this procedure is a convergence of plots which holds 
almost every time we perform the simulation. This turns out to be a consequence 
of the Komlos-Major-Tusnady strong approximation result (Theorem 2 of [23 ): 



sup 
te[o,i] 



S„(t) 



-Z(t) 



o{n- £ ) 



a.s. 



with e 



P 



assuming that -E|£;| p < oo for some p > 3. 



(ii) In our case, the fact that the plots resemble the trajectories of the desired 
limit process suggests that it might be possible to prove some strong approx- 
imation results which would parallel the results presented here. We do not 
investigate this problem here. 



Example 4.3 (Illustration of Theorem 3.7) We consider the coefficients 



j 7 for all j > 1 and Cj = for all j < 0. 



We have two cases: (i) H > 1/a; (ii) H = 1/a. 
To illustrate (i), we let a -1 < 7 < 1 and d n 



,1-7 



i H l / a where 



# = - + 1-7. 

a 



(32) 



The fact that 7 > a 1 guarantees that ( 10 ) holds. Since 



E 



n „ r 

^ ^ - J, 



l- 7 



(1 " 7)- 



the hypotheses of Theorem 3.7 are verified with L2(n) = 1, a 
and 6 = 0. Figure [T] gives an approximation of the sample path of the process 
{■h- a .H,a,o(t)}te[Q,i] obtained for 7 = 0.75; hence a = 4. The picture on the left 
was obtained using a — 1.5 (the limit process is a LFSM with H — 0.92), while 
the picture on the right corresponds to a — 2 (the limit process is the FBM of 
index H = 0.75, multiplied by cl/Ch — 4.28, where Ch is given by p5|.) 

To illustrate (ii), we let 7 > 1 (hence 7 > a -1 ). We consider d n = 1. The 
hypotheses of Theorem |3.7| arc verified since 



n 



d n 



3=1 



Si" 



C(7) 



where C is the Riemann-zeta function. Figure [2] gives an approximation of 
the sample path of the process {aZ(i)} te [ ji obtained for 7 = 4, and hence 
a = ((4) = vr 4 /90 = 1.08 (see p. 807 of [I]). The picture on the left was obtained 



22 



stepfun(time, Z) 



stepfunftime, sum) 




Figure 1: Illustration of Theorem 3.7 for H > a l \ a = 1.5 (left), a — 2 (right) 



stepfunftime. sum) 



stepfunftime, sum) 



00 02 



06 0.8 



0_0 0.2 0.4 0_6 



Figure 2: Illustration of Theorem 3.7 for H — a 1 : a = 1.5 (left), a — 2 (right) 
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using a — 1.5 (the limit is an a-stable Levy motion), while for the picture on 
the right we used a = 2 (the limit is the Brownian motion multiplied by a). 

In both figures, the shift-and-scale of the values /(d n a n ) was done auto- 
matically by the computer, since the convergence holds in the M\ topology. 



Example 4.4 (Illustration of Theorem 3.9) We consider the coefficients: 



hr 1 



if j is even ^ ^ j > 1 and 
it j is odd 



= for all j < 0, 



for some ki,k 2 > 0. We write c 7 



c" where 



o 



and c'j = c" 



if j is even 
if j is odd 



and 



4 = 





k 2 j~ 



if j is even 
if j is odd 



for j > 1 
(33) 



for j < 0. We have two cases: (i) if > 1/a; (ii) if = 1/a. 



To illustrate (i), we let a 1 < 7 < 1 and d n = n 1 7 = n H x l a where if 



is 



given by ( 32 ) . Since 7 > a 



j=X,j even 



[n/2] 
fc=l 



([10]) holds. Note that 
1 



1 



7 



2(1-7) 



j = l,j odd 



Therefore, 



3=1 



J 



=l,j even 



2(1-7) 



3=1 



fcl 



2(1-7) 



and 



1 



2(1-7) 



Hence, a = a! -a" = (k x - fc 2 )/[2(l - 7)]. Note that a = if k x = k 2 . 

Figure[3]gives an approximation for a sample path of the process {A a ^HM,o(t)}te 
obtained for fci = 3, fc 2 = 1, 7 = 0.75 (hence a = 4) and a = 1.5 (left), respec- 
tively a = 2 (right). On the left, the limit process is a LFSM with if = 0.92; on 
the right, the limit is the FBM of index if = 0.75 multiplied by a/C H = 4.28. 
The plot of the values Sj? /(d n a n ) was performed in the range [—4.03,4.83] (left), 
respectively [—1.86,1.75] (right). 

The same values 7 = 0.75 and a = 1.5 were used for Figure [4] but with 
fci = &2 = 1. I n this case, the limit is the zero process since a = 0. In the 
the picture on the left, we imposed the shift-and-scale operation in the interval 
[—0.62,0.52], while in the picture on the right we used the automatic shift- 
and-scale performed by the computer. Therefore, the picture on the right is a 
blow-up of the picture on the left. (Note the small values on the y-axis in the 
picture on the right.) 

To illustrate (ii), we let 7 > 1, d n — 1 and if = 1/a. Then 



1 ™ 



3=1 



E r 

3^1)3 even 



=: a 



and 



1 



Z. c 3 

3=1 



^ E r 

J*>1,3 odd 



=: a 
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stepfunftime, sum) 



stepfun(time. sum) 




T - 



0,0 02 0.4 0.6 0.8 1 .0 0.0 0.2 0.4 0.6 0.8 UB 



Figure 3: Illustration of Theorem 3.9 for k\ = 3, hi = 1, H > a 1 : a = 1.5 
(left), a - 2 (right) 



stepfunftime. sum) 



00 02 



stepfunftime, sum) 



02 0_4 



Figure 4: Illustration of Theorem 3.9 for k\ = k% = 1 (a = 0) and H > a 1 
shift-and-scale in the range [-0.62,0.52] (left), automatic shift-and-scale (right) 
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stepfun(time, sum) 



stepfun(time sum) 




Figure 5: Illustration of Theorem 3.9 for k\ = = 1 5 H = a 1 : a = 1.5 (left), 
a = 2 (right) 

In this case, a — a' — a" — J2j>i c j> but a^O when fci = &2. 

FigurepWives an approximation for a sample path of the process {aZ(t)} te [ ,i] 
obtained for k\ = ki = 1, 7 = 4 and a = 1.5 (left), respectively a = 2 (right). 
In this case, a = J2 j >i(- 1 Vr 4 = -7tt 4 /720 = -0.95 (see p.807 of [lj). The 
picture on the left is an approximation of an a-stable Levy motion, while the 
picture on the right is an approximation of the Brownian motion multiplied by 
a. The plot of the values Su /(d n a n ) was performed in the range [—3.49,4.16] 
(left), respectively [—1.55,1.45] (right). 



Example 4.5 (Illustration of T heorem 3.10) We assume that the coefh 



cients (c,-)j e z are given by (4.4). In order that J2jez\ c j\ < 00 f° r some 
< 5 < a, 5 < 1, we need 7 > max(a , 1). We let 

j£Z j> l,j even i^lj odd 

To illustrate the result, we consider k\ — k-i — \ and 7 = 4 (hence A = 
-0.95). Figure [6] gives and approximation for a sample path of the process 
{AZ(t)} te [ 01 j in the case a — 0.8 (left), respectively a = 1 (right). The plot of 
the values Sj: /a n was performed in the range [—4.38,2.38] (left), respectively 
[-2.32,4.23] (right). 

A Some auxiliary results 

Lemma A.l Let x e D[0, 1] be arbitrary. For any 0<s<u<v<t<l. 
\x(u) - x(v)\ < 2\x(s) - x{t)\ + H(x(s),x(u),x(t)) + H(x(s),x(v),x(t)). 



2G 



stepfun(time. sum) 



stepfun(time. sum) 




Figure 6: Illustration of Theorem 3.10: a = 0.8 (left), a = 1 (right) 



Proof: We consider only the case x(s) < x(t), the case x(t) < x(s) being 
similar. We claim that for any s < u < t, 

\x(u) - x(s)\ < \x(t) - x(s)\ + H(x(s),x(u),x(t)). (34) 

To see this, we consider three cases. If x(u) < x(s), then \x(u) — x(s)\ — 
H(x(s),x(u),x(t)). If x(s) < x(u) < x(t) then \x(u) - x(s)\ < \x(t) - x{s)\. 
Finally, if x(u) > x(t) then \x(u) — x(s)\ < H(x(s),x(u),x(t)) + \x(t) ~ x(s)\. 

The conclusion follows using relation p4| for u and v, and the fact that 
\x(u) — x(v)\ < \x(u) — x(s)\ + \x(v) — x(s)\. □ 

Lemma A. 2 Let x £ D[0, 1] be arbitrary. For < s < t < 1, define 
P= sup H(x(u), x(v), x(w)). 

s<v.<v<w<.t 

Ifr)>2/3 then 

v- P 

where N v (x; [s,t]) denotes the number of n-oscillations of x in the interval [s,i\. 
Proof: Let s < t\ < t2 < <3 < t± < . . . < t2jv-i < t2N < t be such that 

\x(t 2k ) - x(t 2k -i)\ >V for all k = 1, . . . , AT. 
Assume first that xfo) — x(ti) > rj. We claim that: 

^(£3) > x(t 2 ) — (3 and x(ti) — x(t 3 ) > r\. 
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To see this, suppose that x(ts) < x(t 2 ) — (3- Then the distance between x(t 2 ) 
and the interval with endpoints x{t\) and x(t 3 ) is greater than /3, which is a 
contradiction. Hence £(£3) > x(t 2 ) ™ ft- On the other hand, if we assume that 
x(t±) — x(t 3 ) < —77, we obtain that 

x(tt) = x(h) - x(t 2 ) + x(t 2 ) - x(t 3 ) + x(t 3 ) <-ij + I3 + x(t 3 ) < x(t 3 ) - p, 

which means that the distance between x(t 3 ) and the interval with endpoints 
x(t\) and a; (£4) is greater than /3, again a contradiction. 
Repeating this argument, we infer that: 



x{t 2 k) - x(t 2 k-i) > V, f° r a H k = l, 



and 



,N 



,N-l. 



x(t 2k +i) - x{t 2k ) > -P for all k=l, . 
Taking the sum of these inequalities, we conclude that: 

x(t 2N ) - x(h) >Nr)-(N- = N( v -p)+p. 

On the other hand, by Lemma |A.l[ we have: 

\x{t 2N ) - x{h)\ < 2\x(t) - x(s)\ + 2(3. 



Combining (35 1 and (361, we obtain that 

2\x(t) -x(s)\ +/3 



N < 



(35) 
(36) 



v-P 



which is the desired upper bound. □ 
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